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A SIMUTJ^TTON  STUDY  OF  A TWELVE  DEGREE  OF  FREEDOM  SYSTEM 


1 . Introduction 

Boston  College,  In  Its  Investigation  of  the  effects  of  ground  motion 
on  emplaced  missile  systems.  Is  analyzing  various  theoretlrril  models  of 
missile  suspensions  and  their  response  to  motions  received.  The  work 
reported  herein  deals  with  a mathematical  simulation  study  of  a twelve 
degree  of  freedom  system  which  might  be  of  concern. 


1 1 . Mathemat leal  Treatment 

The  Boeing  twelve  degree  of  freedom  model^  for  a missile  suspension 
system  is  a set  of  coupled  ordinary  differential  equations  of  the  second 
order  relating  missile  displacement  to  ground  acceleration.  By  using 
certain  linear  transformations  it  is  possible  to  express  these  equations 
in  terms  of  generalized  coordinates  for  which  the  equations  of  motion 
are  uncoupled.  It  was  desired  to  simulate  this  model  by  numerically 
solving  the  uncoupled  equations  for  applied  accelerations  and  computing 
missile  displacements. 

The  uncoupled  equations  of  motion  are: 

1 f-j  - [rj‘ ' W 

The  matrix  represents  the  ground  accelerations  and  the 

matrix  is  a column  of  displacements  in  generalized  coordinates.  The  other 
matrices  are  functions  of  the  eigenvector  matrix  and  of  the  springs  and 
dashpots  in  the  suspension  system.  The  [^^7  Beatrix  is  diagonal  with  each 
term  equal  to  two  times  the  reciprocal  of  the  decay  time  constant  for  the  mode. 
The  matrix  is  also  diagonal  and  contains  the  elgenfrequencles  of 

the  system  In  the  form  of  the  square  of  the  natural  frequencies  of  each  mode. 
The  [y2  matrix  Is  the  modal  participation  matrix  which  transforms  the 
three  orthogonal  components  of  ground  acceleration  Into  twelve  components 
In  the  generalized  coordinate  system. 

When  this  system  has  been  solved  for  if}  , multiplication  by  the 
eigenvector  matrix  yields  the  relative  displacement  from  the  ground.  Ground 
displacement  may  be  added  to  it  for  absolute  displacement. 
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It  is  believed  that  the  springs  and  dashpots  do  not  remain  constant 
and  It  was  assumed  for  purposes  of  this  simulation  that  the  time  varying 
factors  would  affect  only  the  and  matrices  which  would  remain 

diagonal  but  have  time  varying  elements. 

The  equation  for  the  node  Is: 

-r  ^ ^ +■  ^ where  represents  the 

term.  In  the  matrix  which  Is  the  product  - J- 

Letting  ^ asi)'-  f(t)  we  can 

rewrite  the  equation: 

^ 4-  a,(t)  ^ 

Integrating  once  we  get;  ^ ^ ^ 

assuming  Pit'S- c and  Integrating  a second  time  we  get:  ° . 

/ <■  /tf  * ■>  / ^ 

1-  A*tJ.  ^ tt dr  = J J 

letting  cx.lt)  we  have;^  " * **  ^ 

To  implement  this  in  the  discrete  time  case  we  used  the  sequence  notation 

2 

of  Cuenod  and  DUrllng  : 

A --  [a  , a,,  A,  ^ [al^l  A(i),  a(*X) , a{jAJ.  • • J 

/l8-‘ ,a,4  , ] 

4,  a .f  -a  t a^4-, 

to  obtain  an  approximate  Integral  overtime  using  Trapezoidal  rule: 

■& 

1,1,  •■  • J 

0 

to  approximate  the  second  Integral  by  the  same  method: 

-t-  r 

J falojd^ur  - I,  2,  I ] 

<>  / 

by  applying  this  notation  to  the  Integrated  modal  equation  we  get: 

y^A  A/*[^, '>  V J ■■■]■■  J.  >■  'fi  ■]  -*- 
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culatlons  of  the  twelve  inodes  as  acceleration  and  spring  and  dashpot  values 
change . 
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III.  Validation 

The  driving  force  and  the  damping  coefficient  (a  (t)  ) for  each  mode  was 
made  equal  to  zero  and  an  Initial  displacement  was  given  to  all  the  modes.  The 
time  increment  (h)  was  made  equal  to  .2  seconds  for  observing  the  first  six  modes 
and  .01  seconds  for  the  higher  six  modes.  As  shown  In  Figure  lA  and  Figure  IB 
the  modes  oscillated  at  the  expected  frequencies  without  decaying.  Then  the  a 
(t)'s  were  made  equal  to  the  Boeing-determined  constants  and  the  oscillations 
were  observed  to  be  damped  appropriately  for  each  mode  as  shown  in  Figures  2A 
and  2B. 

It  was  desired  to  run  this  simulation  using  data  for  which  no  vertical 
component  of  acceleration  was  available.  In  order  to  check  the  behavior  of 
the  system  when  subjected  to  vertical  acceleration  a triangular  pulse  of  verti- 
cal acceleration  was  Input  and  missile  motion  plotted  in  Figure  3.  As  seen  in 
this  figure  the  vertical  acceleration  affects  only  the  vertical  component  of 
motion  and  so  It  was  assumed  that  the  rest  of  the  motion  could  be  simulated 
using  only  hte  x and  y components  of  acceleration.  Also  this  data  had  been 
obtained  by  sampling  at  5 samples  per  second,  a rate  too  slow  to  sample  the 
highest  six  modes.  For  this  reason  It  was  decided  to  use  only  the  first  six 
modes.  The  modal  participation  of  the  other  nodes  Is  so  small  as  to  be  negli- 
gible anyway. 

Tlltmeter  data,  scaled  to  represent  acceleration  north  and  west,  was  Input 
to  the  model.  This  data  was  collected  during  a nuclear  explosion  of  magnitude 
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6.3,  1,080  km  away.  Tlip  respon.se  of  the  model  Is  shown  In  Figure  4 which  shows 
missile  motion  with  respect  to  the  silo.  It  shows  that  for  this  event  in  which 
acceleration  reached  a maximum  of  30  micro  g's  peak  to  peak,  the  pitch  term 
reached  a maximum  of  23  arc  seconds  peak  to  peak  and  does  not  reach  a undesir- 
able maximum  whicli  would  occur  about  70  arc  seconds  peak  to  peak. 
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